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ABSTRACT 

The  purpose  of  this  effort  was  to  gain  an  insight  into 
the  Bi-phase  cycle  concept  and  to  develop  an  experimental 
facility  to  test  the  key  component  of  this  power  cycle — the 
Bi-phase  nozzle.   A  system  was  developed  and  checked  out  that 
allows  testing  different  Bi-phase  nozzle  shapes.   The  experi- 
mental system  is  provided  with  sufficient  instrumentation  to 
allow  for  the  deteimination  of  the  nozzle  performance.   In 
addition  to  the  standard  gas  and  liquid  flow  rates,  measure- 
ment provisions  exist  for  survey  of  the  nozzle  exit  plane. 
Measurements  can  be  made  for  the  gas/ liquid  ratio  and 
velocities  at  the  nozzle  exit. 
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I.   INTRODUCTION 

Low  speed  prime  mover  systems  in  marine  applications 
are  attractive  in  that  they  negate  requirements  for  heavy, 
voluminous  and  noisy  reduction  gearing.   Direct  drive  con- 
ventional steam  turbines  and  conventional  electric  drives 
have  been  tested  and  utilized  with  only  limited  success. 
What  is  required  is  a  machine  with  low  specific  speed  that 
is  compatible  at  the  same  time  with  existing  heat  source 
systems.   One  option  that  may  meet  these  requirements  is 
the  Bi-phase  prime  mover.   The  Bi-phase  turbine  is  a 
machine  having  low  specific  speed  and  high  torque.   As 
such,  it  is  ideally  suited  for  marine  propulsion.   The  Bi- 
phase turbine  can  drive  the  propellers  directly  without 
gear  reduction.   This  feature  leads  to  improvements  in 
weight,  volume  and  noise  reduction. 

In  addition,  certain  aspects  of  the  Bi-phase  concept 
may  allow  increased  power  production  efficiencies.   It  may 
be  possible  to  achieve  near  constant  temperature  expansion 
in  the  two-phase  nozzle  thus  realizing  a  degree  of 
"Carnotization"  of  the  basic  Rankine  cycle.   A  detailed 
description  of  the  Bi-phase  power  plant  is  given  in  Section  III. 

A  key  element  in  the  development  of  the  Bi-phase  prime 
mover  is  the  Bi-phase  nozzle.  A  Bi-phase  nozzle  is  a  flow 
accelerator  in  which  fluid  particles  or  droplets  are 
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accelerated  by  a  vapor.   The  liquid  and  vapor  may  be 
chemically  identical  or  dissimilar.   The  expanding  vapor 
imparts  kinetic  energy  to  the  liquid  droplets.   During 
this  expansion,  the  droplets  may  also  transfer  heat  to  or 
from  the  vapor,  resulting  in  a  near  isothermal  expansion. 
In  the  nozzle,  most  of  the  available  thermal  energy  of  the 
vapor  is  converted  to  kinetic  energy  of  the  liquid  droplets. 
The  nozzle  is  converging-diverging,  a  requirement  resulting 
from  the  low  sonic  velocity  in  two-phase  mixtures.   Computer 
codes  [Ref.  1]  exist  that  yield  nozzle  exit  velocities  to 
reasonable  accuracies  provided  fluid  properties  are  known, 
liquid  particle  size  can  be  controlled,  and  one  dimension- 
ally  in  the  flow  field  is  assumed. 

Tangren,  Dodge,  and  Siefert  derived  the  equations  for 
isentropic  flow  of  an  immiscible  two-component  mixture 
with  equal  temperatures  of  the  liquid  and  the  gas.   The 
first  extensive  experiments  on  two-phase  nozzles  were 
performed  by  Elliott  [Ref.  1]. 

The  purpose  of  this  effort  is  to  develop  an  experimental 
facility  that  will  allow  evaluation  of  two-phase  nozzle 
concepts  and  associated  liquid-vapor  interaction  phenomenon. 
The  ultimate  goal  is  to  provide  experimental  capability  to 
validate  analytic  models  of  Bi-phase  nozzles.   To  this  end 
an  experimental  system  was  designed  that  includes  a  high 
degree  of  flexibility  to  allow  for  rapid  flow  passage 


geometry  variation.  The  instrumentation  installed  has  the 
capability  to  sample  liquid  content  at  different  points  in 
the  nozzle  exit  plane. 
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II.   BACKGROUND 

Turbine  type  prime  movers,  whether  based  on  the  Rankine 
or  Brayton  cycles,  have  many  advantages.   Among  these  are: 
high  torque  at  low  speed,  low  vibration  levels  as  compared 
to  internal  combustion  energy  converters,  and  in  some  cases 
very  high  power  densities.   Practical  considerations  result 
in  prime  mover  power  outputs  which  are  at  relatively  high 
rotational  speeds. 

Power  demand  devices  such  as  propellers  and  pump  jets 
have  rotational  speed  requirements  set  by  considerations 
considerably  different  than  those  of  the  prime  mover.   Un- 
fortunately the  rotational  speeds  of  prime  mover  and  power 
user  seldom  match.   As  a  consequence,  a  power  conditioning 
device,  a  gear  box  or  reduction  gearing  is  invariably 
mandated.   In  certain  installations  the  reduction  gearing 
may  weigh  more  than  the  prime  mover. 

From  a  purely  thermodynamic  consideration  one  would 
want  the  maximum  cycle  temperature  to  be  as  high  as  possi- 
ble.  However,  where  the  working  substance  must  be 
contained  by  physical  structures  such  as  combustors, 
boilers,  nozzles,  and  turbine  blades,  the  cycle  temperature 
limit  is  set  by  material  considerations.   Furthermore, 
secondary  aspects  of  a  prime  mover,  i.e.,  safety,  may 
further  restrict  the  range  of  cycle  parameters.   A  case  in 
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point  is  the  Navy  nuclear  power  plant  that  does  not  admit 
superheat.   The  turbine  inlet  conditions  are  those  of 
saturated  vapor.   This  constrains  the  maximum  cycle  effi- 
ciency to  comparatively  low  values. 

A  relatively  new  prime  mover  concept  that  may  have  the 
capability  to  circumvent  some  of  the  above  identified  diffi- 
culties is  the  Bi-phase  prime  mover.   The  technology  base 
for  the  Bi-phase  concept  is  founded  in  earlier  work  at  the 
Jet  Propulsion  Laboratory  on  liquid  metal  power  generating 
concepts. 

The  baseline  Bi-phase  concept  is  still  the  Rankine 
cycle  (Fig.  1) .   The  differences  are  subtle  but  of  extreme 
importance.   The  working  fluid  is  vaporized  in  the  heat 
addition  component  of  the  cycle.   This  may  be  a  conventional 
boiler  generating  high  temperature  and  pressure  steam. 
Normally  the  working  fluid  would  then  enter  a  nozzle/turbine 
system  and  energy  would  be  converted  to  a  mechanical  form. 
In  the  Bi-phase  system,  however,  a  second  fluid  is  intro- 
duced at  boiler  temperature  and  pressure  just  prior  to  the 
point  where  the  steam  enters  the  nozzles.   This  second 
fluid  may  have  a  vaporization  temperature  higher  than  water 
and  thus  will  remain  in  the  liquid  state  throughout  the 
cycle.   In  the  nozzle  proper  the  vapor  constituent  acceler- 
ates and  in  so  doing  drags  the  droplets  of  the  liquid  phase 
to  higher  velocity  (Fig.  2) .   This  energy  exchange  of  enthalpy 
drop  of  the  vapor  to  kinetic  energy  of  the  liquid  is  the 

major  departure  from  the  conventional  steam  power  plant. 

12 


BI-PHASE  CYCLE 


400  r- 


ENTROPY  (kj/kg) 


Fig.  1,   Temperature-Entropy  Diagram 
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Fig.  2.   Liquid  Entrainment  Photograph 
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Downstream  of  the  nozzle  the  high  kinetic  energy  liquid 
and  the  expanded  vapor  are  separated.   The  vapor  is  condensed 
and  returned  to  the  boiler.   The  liquid  is  "parked"  on  a 
rotating  wheel  and  energy  is  extracted  by  an  impulse  or  pitot 
type  turbine.   After  interaction  with  the  hydraulic  turbine, 
the  fluid  is  returned  for  injection  at  the  Bi-phase  nozzle 
mixing  section. 

It  is  the  Bi-phase  aspect  of  the  system  that  is  at  the 
base  of  some  of  the  potential  advantages  of  the  concept, 
which  are: 

1.  The  energy  imparted  to  the  liquid  phase  is  in  a 
kinetic  form  of  a  high  density  mass.   Hence,  the 
liquid  velocities  involved  are  relatively  low  as 
compared  to  velocities  of  the  vapor.   Thus  the 
output  of  the  impulse  hydraulic  turbine  will  be  at 
high  torque  and  low  rotational  speeds — more  in 
keeping  with  rotational  speeds  of  driver  elements. 

2 .  The  energy  extraction  from  the  steam  and  the  energy 
addition  to  the  liquid  takes  place  in  an  inherently 
small  device — the  Bi-phase  nozzle.   Thus,  the  high 
temperature,  pressure,  and  velocities  are  confined 
to  a  high  energy  density  component. 

3.  Perhaps  the  most  important  aspect  is  the  phenomenon 
that  takes  place  in  the  nozzle  proper  between  the 
liquid  and  vapor  and  how  this  interaction  impacts 
the  overall  cycle  performance.   While  it  is  normal 
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practice  to  assume  "isentropic"  expansion  in  a 
nozzle,  it  is  not  likely  that  this  is  what  happens 
in  a  Bi-phase  expansion.   In  all  probability  thermal 
energy  is  exchanged  between  the  liquid  and  vapor 
constituents.   This,  in  fact,  may  be  the  most 
attractive  aspect  of  the  Bi-phase  concept.   If 
thermal  energy  can  be  added  to  the  vapor  phase 
from  the  liquid,  then  the  expansion  can  be  made  to 
approach  an  "isothermal"  process.   The  net  result 
is  a  degree  of  "Carnotization"  of  the  basic  Rankine 
cycle  that  can  only  lead  to  increased  overall  cycle 
efficiencies. 
Application  of  two-phase  turbines  to  Navy  propulsion  has 

been  examined  in  two  design  studies.  References  2  and  3. 

In  both  cases,  the  following  potential  advantages  were  found 

when  the  Bi-phase  turbine  was  compared  to  advanced  steam 

turbine  systems. 

1.  High  efficiency — Full  load  output  power  performance 
gains  ranging  from  2  0  to  50  percent  were  found. 

2.  Direct  drive — Direct  drive  at  speeds  ranging  from 
90-4500  rpm  was  found  possible  with  a  single-stage 
turbine. 

3.  Reduced  volume — Volume  reductions  of  30  percent 
were  estimated. 

4.  High  part-load  efficiency — Variable  mass  ratio 
enabled  part-load  (cruise)  efficiency  gains  of  as 
much  as  100  percent. 
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An  independent  Navy  study  of  a  large  Bi-phase  marine 
turbine  arrived  at  these  same  conclusions.   These  potential 
advantages  lead  to  many  possible  applications  for  marine 
use.   Most  prominent  among  these  advantages  are  those 
applications  where  space  and  weight  are  the  most  important 
considerations,  for  example,  torpedoes,  submarines,  SES, 
SWATH,  etc. 

It  is  thus  apparent  that  the  key  element  of  a  Bi-phase 
power  system  is  the  Bi-phase  nozzle.   It  is  this  component 
that  is  the  subject  of  this  thesis.   Specifically,  the  goal 
is  to  develop  a  research  capability  to  study  the  flowfield 
and  energy  interaction  in  a  Bi-phase  nozzle. 
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III.   BI-PHASE  CYCLE  THERMODYNAMICS 

The  Bi-phase  cycle  is  basically  a  variation  of  the  con- 
ventional Rankine  steam  plant.   In  reality,  the  "Bi-phase 
aspect"  should  be  a  inulticomponent  cycle  as  there  is  no 
requirement  that  the  components  (vapor  and  liquid)  of  the 
working  substance  be  chemically  identical. 

In  a  conventional  Rankine  cycle  (Fig.  1)  the  liquid 
phase  is  introduced  to  the  heat  additive  component  of  the 
system  at  the  maximum  cycle  pressure.   Energy  is  added  and 
the  substance  is  vaporized.   This  is  a  constant  pressure 
energy  addition  process  and  follows  a  constant  pressure  line 
on  the  T-S  coordinate  system.   The  maximum  cycle  temperature 
may  be  well  into  the  superheat  region.   Expansion,  often 
idealized  as  isentropic,  follows.   During  this  pressure  drop 
from  boiler  to  condenser  pressures,  work  is  extracted  by  a 
turbine  via  a  volume  expansion.   Constant  pressure  heat 
rejection  returns  the  working  substance  to  a  liquid  phase 
which  then  is  pumped  up  to  boiler  pressure  and  reintroduced 
to  the  heat  addition  component. 

A  common  way  of  depicting  the  maximum  cycle  efficiency 
for  any  cycle  is  to  apply  the  Carnot  efficiency  criteria. 
This  Carnot  efficiency  depends  only  on  temperatures  of  heat 
addition  and  heat  rejection. 
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,    Tc 
Hn  =  1  -  T^ 

Where  n   =  Carnot  efficiency 

T^  =  Condenser  temperature  ('-'K) 

T  =  Boiler  temperature  (°K) 
o 

Thus,  a  cycle  operating  between  source-sink  temperatures  of 
232*^C  and  64*^0  respectively  has  a  maximum  energy  conversion 
of  33%.   The  efficiency  for  a  Rankine  cycle  with  the  same 
temperature  limits  and  saturated  vapor  turbine  inlet  condi- 
tions yields  an  ideal  efficiency  of  29%.   This  efficiency 
is  based  strictly  on  thermodynamic  considerations  and  does 
not  consider  irreversibilities  such  as  friction  losses, 
temperature  drops  in  heat  exchangers,  and  hydrodynamic  losses 
in  the  cycle  mechanical  components.   Furthermore,  from 
practical  considerations,  the  last  stages  of  the  expander 
see  a  very  wet  working  fluid  =  75%  quality.   This  results 
in  erosion  and  impact  losses  leading  to  further  decrease 
in  the  power  generator  effectiveness. 

It  is  the  Bi-phase  cycle  variation  that  has  the  potential 
to  circumvent  the  above  identified  difficulties  of  the 
Rankine  plant. 

The  Bi-phase  cycle  adhieves  improved  efficiencies  because 
the  process  involved  allows  the  development  of  a  near  Carnot 
cycle.   Thus,  the  Bi-phase  concept  achieves  a  degree  of 
Carnotization  of  the  Rankine  cycle.   This  Carnotization 
takes  place  in  the  Bi-phase  nozzle  process.   While  in  a 
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Rankine  plant  the  expansion  process  is  idealized  as  isen- 
tropic;  in  the  Bi-phase  system,  it  can  be  caused  to  approach 
an  isothermal  process. 

There  are  basically  two  Bi-phase  techniques  whereby  a 
Rankine  cycle  can  be  caused  to  approach  the  Carnot  cycle. 

A.   TWO  COMPONENT  SYSTEM 

One  method  to  Carnotize  a  steam  cycle  is  to  add  a  second, 
chemically  different  component  that  remains  a  liquid  through- 
out the  cycle  process.   This  component  provides  a  source  of 
heat  for  the  expanding  vapor  producing  a  near- isothermal 
expansion.   In  Fig.  1  it  can  be  seen  that  the  temperature  at 
the  end  of  the  nozzle  expansion  is  218  C  compared  to  64  C 
which  would  occur  in  the  Rankine  cycle.   This  results  in 
a  cycle  efficiency  of  31%  versus  29%  for  the  Rankine  cycle. 
This  cycle  was  invented  by  D.  G.  Elliott  of  Cal  Tech  for 
liquid  metal  MHD.   [Ref.  1] 

The  two-component  cycle  uses  a  low  vapor  pressure  liquid 
and  a  high  vapor  pressure  liquid  having  different  chemical 
makeup.   Some  fluid  combinations  which  have  been  considered 
are  steam-Krytox,  steam-Caloria,  steam-lead  bismuth  eutectic, 
and  Dow-therm-Therminol.   The  basic  advantages  of  a  two- 
component  engine  are  low  rpm,  high  efficiency  with  no  steam 
extraction  and  a  more  compact  turbine.   Disadvantages 
include  larger  heat  exchangers  and  the  complexity  of  handling 
two  fluids.   Figs.  3  and  4  depict  two  possible  methods  of 
utilizing  the  two-component  approach. 
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Fig.  3.   U-Tube  Turbine  Engine 
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Fig,    4.      Impulse   Turbine    Engine 
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D.   SINGLE  COMPONENT  SYSTEM 

In  this  cycle  the  expansion  process  is  ideally  paralled 
to  the  saturated  liquid  line  that  outlines  the  left  hand 
side  of  the  enthalpy  dome.   In  order  to  achieve  this  type 
of  Carnotizing,  heat  has  to  be  removed  during  expansion. 
While  it  may  be  possible  to  do  that  through  the  walls  of 
the  nozzle,  a  more  conventional  way  is  to  extract  steam  in 
steps.   That  means  that  a  scheme  of  pressure  staging  is 
required,  where  after  each  expansion  stage  steam  is  ex- 
tracted and  used  to  preheat  feed  water  by  condensing  all  of 
the  extracted  steam. 

The  single  component  system  has  the  advantage  of  elimin- 
ating oil  or  heat  transfer  fluid  as  the  second  component. 
The  two-phase  mixture  in  this  case  consists  of  steam  and 
water.   Typical  volume  ratios  of  water  to  steam  are  10:1 
or  less  so  this  cycle  produces  somewhat  higher  values  of 
turbine  speed. 

In  order  to  obtain  high  efficiency  several  expansion 
stages  are  used.   The  five  stage  system  resulting  in  the 
temperature-entropy  variation  of  Fig.  5  is  shown  schematic- 
ally in  Fig.  6.   Steam  from  the  heat  exchanger  flows  to  the 
first  stage  nozzle.   A  portion  of  the  steam  is  also  used 
to  heat  the  water  flowing  to  the  first  stage  nozzle.   The 
mixture  of  primary  steam  and  water  is  expanded  in  the 
nozzle  from  an  inlet  vapor  quality  of  15%  to  an  exit 
quality  of  18%.   The  two-phase  mixture  drives  a  rotary 
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separator  turbine.   Part  of  the  separated  steam  is  extracted 
and  flows  to  a  feed  water  heater.   The  remainder  of  the 
steam  is  remixed  with  the  water  exiting  the  first  stage 
turbine.   This  mixture  is  expanded  in  the  second  stage 
nozzle  from  an  inlet  vapor  quality  of  12%  to  an  exit 
quality  of  18%. 

The  two-phase  mixture  drives  the  second  stage  rotary 
separator  turbine,  where  additional  steam  extraction  occurs. 
In  the  final  stage ,  steam  exhausting  from  the  turbine  flows 
to  the  condenser.   The  condensate  is  remixed  with  the 
liquid  exiting  the  last  stage  turbine.   The  feed  water 
flows  through  several  stages  of  a  feed  water  heater  where 
the  extraction  steam  is  condensed  on  the  feed  water. 

Five  stages  are  shown  in  Fig.  6.  The  actual  number  of 
stages  will  be  determined  by  a  tradeoff  between  efficiency 
and  cost  and  complexity. 
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IV.   BI-PHASE  NOZZLE  THEORY 

The  Bi-phase  nozzle  is  the  key  element  in  the  Bi-phase 
engine  concept.  While  initially  it  may  appear  that  the 
design  of  the  nozzle  profile  is  quite  straight  forward, 
closer  examination  will  reveal  that  the  question  is  con- 
siderably more  complex.  The  complexity  steams  from  the 
nature  of  the  working  fluid  itself.  It  is  a  mixture  of 
liquid  particles  and  a  compressible  gas  or  vapor. 

The  complicating  factors  can  be  grouped  into  three 
categories : 

1.  The  working  fluid  has  two  constituents  significantly 
different  in  density.   Thus,  there  is  an  aerodynamic 
interaction  between  the  components  (liquid  droplets 
are  accelerated  by  the  vapor) . 

2.  The  shape  and  size  of  the  liquid  particles  is 
governed  by  the  forces  on  them.   In  this  particular 
case  the  surface  tension  and  inertia  forces  are  of 
key  interest. 

3.  Thermal  interaction  between  the  working  fluid 
components  can  range  from  none  at  all  to  phase 
change  of  the  liquid  to  its  vapor  foinn. 

Another  factor  of  interest  in  Bi-phase  flow  is  the 
question  of  the  velocity  of  propagation  of  a  disturbance, 
i.e.,  velocity  of  sound.   While  the  velocity  in  a  particular 
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vapor  and  in  a  particular  liquid  are  well  established,  the 
situation  in  a  mixture  of  vapor  and  liquid  is  considerably 
different. 

Consider  the  following  discussion  from  Olson  [Ref.  4]. 
The  acoustic  velocity  in  any  medium  is  c  =  /K/p,  where  K  is 
the  modulus  of  elasticity  for  the  medium  and  o  is  the 
density.   For  gas-liquid  mixtures  the  acoustic  velocity 
becomes  less  than  that  for  either  the  liquid  or  gas  alone. 
For  a  liquid  with  a  small  concentration  of  gas  nuclei,  the 
elastic  modulus  of  the  mixture  is  reduced,  with  no  appre- 
ciable reduction  in  density,  and  thus  the  acoustic  velocity 
is  reduced.   For  a  gas  with  minute  liquid  droplets,  the 
density  of  the  mixture  is  increased  with  no  appreciable 
change  in  elastic  modulus,  and  again  the  acoustic  velocity 
for  the  mixture  is  reduced.   The  velocity  of  sound  in  a 
14  percent  water  and  86  percent  air  mixture  is  about  30  m/s 
(depending  on  the  pressure,  temperature,  and  impressed 
frequency)  as  compared  with  about  1460  m/s  for  water  and 
about  335  m/s  for  air. 

The  flow  phenomenon  of  a  Bi-phase  mixture  has  been 
analyzed  in  Ref.  1.   It  is  repeated  as  follows.   The 
problem  is  illustrated  in  Fig.  7.   A  spatially  uniform  two- 
component  mixture  of  liquid  drops  and  gas  enters  a  nozzle 
at  high  pressure  and  low  velocity  and  expands  to  low 
pressure  and  high  velocity.   The  objective  of  the  analysis 
is  to  determine,  for  a  specified  pressure,  the  drop  diameter 
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Fig.  7.   Two  Phase  Nozzle  Flow  and  Nomenclature 
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D  and  the  temperatures  T   and  T^ ,    velocities  V  and  V_  and 

g      L  g      L 

flow  rates  m  and  m..  of  the  gas  and  liquid  phases,  respec- 
tively, at  each  station  in  the  nozzle  given  the  initial 
values  of  D,  T^,  T^.  ,  V^,  V^.  ,  the  total  flow  rate,  and  the 

g      Li      g      Li 

properties  of  the  fluids. 

The  five  relations  employed  to  compute   the  five  unknowns 
D,  T  ,  T- ,  V  ,  and  V- ,  are  (1)  the  momentum  equation  for  the 
mixture,  (2)  the  energy  equation  for  the  mixture,  (3)  the 
drop  drag  equation,  (4)  the  drop  heat  transfer  equation,  and 
(5)  the  drop  breakup  criterion.   Solubility  and  vapor 
pressure  relations  provide  the  flow  rate  ratio  m  /m  . 

A.   ASSUMPTIONS 

The  assumptions  employed  in  the  two-component  analysis 
are  as  follows: 

1.  The  liquid  is  uniformly  dispersed  as  spherical  drops 
all  of  the  same  diameter. 

2.  The  drops  break  up  to  limit  the  Weber  number  to  6. 

3.  There  are  no  external  forces  acting  on  the  two-phase 
mixture  other  than  pressure  and  wall  shear,  and 
there  is  no  heat  transfer  to  or  from  the  mixture. 

4.  The  flow  is  one-dimensional. 

5.  The  drops  are  large  enough  for  the  surface  curvature 
to  have  negligible  effect  on  the  vapor  pressure  of 
the  liquid  and  for  the  surface  energy  to  be  negligible 

6.  The  drops  are  isothermal. 

7.  The  gas  mixture  obeys  the  additive-pressure  law. 
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8.  The  partial  pressure  of  the  predominantly  liquid 
component  in  the  liquid  is  given  by  Raoult's  Law. 

9.  The  concentration  of  the  predominantly  gaseous 
component  in  the  liquid  is  given  by  Henry's  Law. 

10.   The  volume  of  the  liquid  solution  is  equal  to  the 

sum  of  the  volumes  of  the  pure  liquids. 

Assumption  1  restricts  the  analysis  to  nozzles  having 

spatially  uniform  injection  of  the  liquid  into  the  gas  and 

operating  at  gas-to-liquid  volume  ratios  greater  than 

unity.   Assumption  2,    the  drop  breakup  criterion,  states 

that  drop  diameter  is  limited  to  a  value  D  for  which  W 

W  =  p  V^  D/2a  =  6.   Thus 
e   "^g  g  ^ 

D    =  ^^^  (1) 

°max      2  (1) 

g  g 

where  p   is  the  gas  density,  V   is  the  slip  velocity  V  -  V^ , 
g        ^         "^    s  ^        "^   g   L 

and  a  is  the  liquid  surface  tension.   The  form  of  Eq.  (1)  is 
physically  reasonable  in  that  the  Weber  number  is  proportional 
to  the  ratio  of  stagnation  pressure  p  Vg/2  to  surface 
tension  pressure  4a/D.   Hence,  a  drop  would  be  expected  to 
flatten  and  break  up  at  a  sufficiently  high  value  of  We. 
This  has  been  verified  experimentally  and  the  critical  Weber 
number  found  to  be  6,  within  a  factor  of  about  two.   An 
additional  restriction  is  that  for  actual  breakup  to  occur, 
the  time  spent  at  a  Weber  number  exceeding  6  must  be  longer 
than  the  natural  period  of  oscillation  of  the  drop. 
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3   h 

tt(Pt.D  /a)  /4,  where  p^.  is  the  density  of  the  liquid.   This 

Li  Li 

requirement  is  met  only  in  two-phase  nozzles  longer  than 
about  10  in.  and  Assumption  2  may  cause  the  analysis  to 
overestimate  the  exit  velocity  by  increasing  amounts  as 
the  nozzle  length  decreases  below  10  in. 

Assumption  3  excludes  magnetohydrodynamic  and  mechanical 
body  forces.   The  exclusion  of  wall  heat  transfer  is  correct 
for  the  insulated  nozzles  of  interest  for  power  systems. 
In  addition,  the  relatively  high  velocity  results  in  short 
residence  lines  in  the  nozzle  proper. 

Assinnption  4  is  closely  met  in  practical  nozzles  since 
good  performance  requires  small  wall  angles,  large  throat 
radius  of  curvature,  and  uniformly  distributed  injection 
of  the  fluids  at  the  nozzle  entrance. 

Assumption  5  is  valid  for  the  drop  sizes  of  0.001  to 
0.010  in.  produced  by  the  Eq.  (1)  breakup  criterion. 
Assumption  6  is  valid  because  of  the  rapid  internal  circula- 
tion in  drops.   Assumption  7  introduces  negligible  error 
in  most  cases  of  practical  interest  since  the  vapor  pressure 
of  the  liquid  is  small  and  needs  only  to  be  evaluated 
approximately . 

Assumptions  8,  9,  and  10  are  either  valid,  or  cause 
little  error,  for  fluids  of  low  miscibility,  which  are  the 
fluids  of  interest. 
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B.   DERIVATION  OF  EQUATIONS  FOR  FREE-STREAM  FLOW 

1.  Continuity 

Referring  to  Fig.  1,  the  nozzle  flow  area  A  is  equal 
to  the  gas  flow  area  m  /o  V  plus  the  liquid  flow  area 

g  '  g  g 

m^/p^V^.   Thus, 

^  (^Pg^g    Pl^l 

where  r  is  the  mass  mixture  ratio  nL/m  . 

2 .  Momentum 

By  Assumption  3,    the  only  force  acting  on  the  free- 
stream  flow  is  that  due  to  the  pressure  gradient.   If  M  is 
the  momentum  flux  at  flow  area  A,  the  change  in  momentum 
flux  across  pressure  increment  dp  is 

dM  =  -Adp  (3) 

The  momentum  flux  can  be  written  as  the  sum  of  the 
momentum  fluxes  of  the  gas  and  liquid.   Thus, 

M  =  m  V   +  m^V,  (4) 

g  g    1  1 

If  the  flow  were  allowed  to  continue  at  constant 
pressure,  V  and  V  would  become  equal  to  each  other  at  the 
mass-weighted  mean  velocity  V.   Since  for  this  process, 
dM  =  0,  the  value  of  V  is  given  by 


(m^  +  m. )V  =  m  V^  +  m^ V^  (5) 

g     i.      g  g     11 
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or 

V  +  rV. 


V  =  -^^r—r — -  (6' 

1  +  r 


Thus,  the  momentum  flux  can  be  written 


M  =  (m   +  m,  )V  (7) 

g     1 

Since  m  +  m^  is  constant,  the  change  in  momentum  flux  is 

g        Li 

dM  =  (m  +  m, )dV  (8) 

g    1 

Substituting  Eqs.  (8)  and  (2)  into  Eq.  (3) ,  dV  can  be  written 
The  slip  ratio  is  defined  as 


s  =  V  /V  =  (V   -  V, )/V  (10) 

g        g     X 

This  equation  can  be  combined  with  Eq.  (6)  to  give  V  and 
V_  in  terms  of  V: 

\  =  (i  -^  T^  1  V  =  av  (11) 


bV  (12) 

The  gas  density  can  be  expressed  as 


p   =  W  p/RT  (13) 

^g    g^^   g 

where  W  is  the  effective  molecular  weight  of  the  gas  mixture 
g 

and  R  is  the  universal  gas  constant.   Eq.  (13)  is  the 
definition  of  the  effective  molecular  weight  W  ,  which  is  the 

g 
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quantity  that  gives  the  actual  gas  density  when  substituted 
in  Eq.  (13) . 

Substituting  Eqs.  (11) -(13)  into  Eq.  (9),  the 
differential  momentum  equation  is 

The  quantities  a  and  b  are  slowly  varying  because  s  is 
typically  only  0.1  to  0.3  and  slowly  varying.   The  quanti- 
ties r,  T  /  W  and  p^  are  also  slowly  varying.   Integrating 
Eq.  (14)  over  a  pressure  increment  Ap,  for  which  a,  b,  r, 

T  ,  W  ,  and  p_  are  constant  to  within  the  desired  accuracy, 

2 

the  change  m  V  is 


(15) 


All  quantities  other  than  pressure  can  be  taken  outside  the 
integral  and  evaluated  at  their  mean  values  (denoted  by 
subscript  m)  corresponding  to  the  mid-interval  pressure  p. 
Thus, 


AV2  =  - 


dp  (16) 
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Performing  the  integrations. 


RT 


AV2  =  - 


'm 


1  +  r. 


'm 


a  W 
m  g 


log 


P  +  Ap/2 


r  Ap 


e  p  -  Ap/2    b^pl 


(17) 


m 


m  m 


Equation  (17)  is  the  final  form  of  the  momentum  equation. 
3.   Energy 

The  enthalpy  change  of  the  mixture  between  state  1 
(the  beginning  of  pressure  interval  Ap)  and  state  2  (the 
end  of  the  interval)  can  be  evaluated  in  two  steps: 


(1)   phase  change  at  p,  ,  T   ,  T,.   and  (2)  change  to  P^ ,  T 

"l  1 

fixed  composition. 

The  enthalpy  change  for  step  1  is 


T  2 /  ^t  fixed  composition. 


AH^  = 


+ 


amount  of  Ai 
vaporized   J 


amount  of  B 
vaporized 


amount  of  A 
and  B 
vaporized 


X 


X 


X 


enthalpy  required 

to  vaporize  and 

heat  unit  mass  of 

A  from  Tt  to  T  ^ 

L     gL  ' 

enthalpy  required 
to  vaporize  and 
heat  unit  mass  of 
B  from  V-   to  T 


kinetic  energy 
required  to 
accelerate  unit 
mass  from  V   to 
V  ^1 

^1 
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or 


AH,  =  (m^    -  m^   )  [L    +  c^   (T    -  T,  )  ] 

g-,    go   1    gi   1    i 


'2      ^2     "      ^1 

-  A   )  (V2   -  vf 


+  (iTi^   -  A   )  (V2   -  vf  )/2  (18) 

g^   gT   g-    i- 


where  L  and  c  are  latent  heat  and  specific  heat,  respectively, 
Introducing  more  compact  notation. 


AH,  =  Am   (L    +  c    5,T) 
1     a   a^     a   1 

g   1    g^ 

Am  5,v2 

+  Amj^  {L^      +  c^   5^T)  +  — 2^i (19) 

2    1      g^ 


The  enthalpy  change  for  step  2  is  evaluated  from  the 
temperature,  pressure,  and  velocity  changes,  with  properties 
evaluated  at  mean  T  and  p  for  the  interval. 

AH^  =  m^  [c^  (T    -  T^  )  +  y  {vl      "  v2  ) ] 
2    ^2   ^m  ^2    ^1    2   g2    g^ 


+  m,  [c,  (T,   -  T,  )  +   ^     ^ 


^2   ^m   ^2     h  Pi 


m 


+  ^  (V2   -  vf  )  ]  (20) 

"^    -^2     ^1 
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AV 


=  m   (c   AT   + 


+  m   (c   AT   +  ^^£-  +  -y 

^2V  m   -^    pi      ^ 
m 


(21) 


3y  Assumption  3,  no  work  is  done  by  the  free-stream  flow  and 
no  heat  is  transferred  to  it.   Hence, 


AH   +  AH2  =  0 


(22) 


Substituting  Eqs.  (19)  and  (21)  into  Eq.  (22)  and 

solving  for  AT   gives  the  energy  equation  for  the  mixture; 
y 


AT   = : 


Av2 


'm 


\  m      pi 

^   m 


Am  6,v2    ^*a 


2ih 


m 


(L    +  c^    5,T) 

^  ^1 


Am, 


+  -.— 2  (L,   +  c,   6^T) 
m      b,     b^   l 
gn     1     ^1 


(23) 


4.   Drag 

Although  no  force  other  than  pressure  acts  on  the 
free-stream  flow  as  a  whole,  a  drag  force  exists  between 
the  phases.   Hence,  a  second  momentum  equation  must  be 
written  using  as  the  control  volume  the  boundary  between 
the  phases. 

The  two  forces  acting  on  each  liquid  drop  are  the 
buoyancy  due  to  the  pressure  gradient  and  the  drag  due  to 
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the  relative  gas  velocity.   The  sum  of  these  is  equal  to 
the  mass  times  the  acceleration  of  the  drop.   Thus,  for  a 
single  drop. 


dynamic  pressure  of 
relative  gas  flow 


X 


drag  co- 
efficient 


X 


frontal  area 
of  drop 


volume 
of  drop 


X 


pressure 
gradient 


mass 
of   drop 


X 


accelera- 
tion of   drop 


or 


(l^gl^2l^2)    SID2    _^dH 


^.d3       \/^      ^\\ 
-6-   %'  Vl   ^J 


(24) 


The  absolute  value  sign  in  the  first  term  makes  the 
drag  force  positive  when  V  >  V_  and  negative  when  V  <  V  • 
Solving  Eq.  (24)  for  dV  , 


dv,  = 


Sp^lslsV^C^dx 

g'  ' P 

4p^V^D 


_  _^ 


P^V, 


(25) 


Differentiating  Eq.  (12) ,  dV  can  also  be  expressed  in  terms 
of  s,  r,    and  V.   Thus, 


—   —  f    sdr       ds  1 

dv^  =  bdv  +  V  [jTir-^yz  "  FTT  ! 


(26) 


Solving  for  ds. 


T-/TI  \  ^^             J  (1  +  r)  dV, 

J    b(l  +  r)dV  ,  sdr  1 

ds  =  +  ^i ; -  

V  ^  ^  ^  V 


(27) 
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Substituting  dV   from  Eq.  (25) ,  noting  that  dV  =  dv2/2V, 
using  Eq.  (12) ,  and  writing  for  a  finite  increment,  results  in. 


As  = 


b  (1  +  r  ) Av2    (1  +  r  )Ap    s  Ar 
m m       nL_  +    ^ 


2v2 
m 


b  p,  v2 
m'^1   m 
m 


1  +  r 


m 


3p^  |s^|s  C^  (1  +  r^) Ax 
g  '  m '  m  D       m 
^m m 

mi 
m 


(28) 


This  is  the  drag  equation  employed  when  x  is  specified 
as  a  function  of  p. 

Solving  Eq.  (28)  for  Ax  yields  the  required  alterna- 
tive equation: 


Ax  = 


4D 


3p^   s^  s  C^  V2 
g  I  m  I  m  D 
^m'   '    m 


Ap  + 


b^p,  AV2 
ml 
m 


1 


m'^l   m  /  s  Ar 
+   ,  .^   (  ^ -  Asjl   (29) 


1  +  r^  V 1  +  r 
m  ^ 


m 


5 .   Heat  Transfer 

Although  no  heat  is  transferred  to  the  mixture  as 
a  whole,  heat  transfer  exists  between  the  phases.   Hence, 
a  second  energy  equation  must  be  written  using  as  the 
control  volume  the  boundary  between  the  phases. 

The  work  dW  done  on  the  liquid  is  that  due  to  drag 
by  the  gas.   (Only  work  done  by  shear  or  shaft  forces  is 
included  in  dW  when  writing  the  First  Law  for  a  control 
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volume) .   Multiplying  Eq.  (24)  by  the  number  flow  rate  of 

3 

i_/TrD  p,,  the  drag  force  F, 


3 
drops  N  =  6m_/TrD  p,  ,  the  drag  force  F,  on  that  quantity  of 


liquid  is 


•  m^  ,         dV^ 

F.  =  ^P^|V  iVC^TTD^  =  _i  §£  +  m^V^  -3-^  (30) 

d   8*^g'  g'  g  D       ^  dx    1  1  dx 


The  work  done  on  the  liquid  is 

-dW  =  F,dx  =  m  (^  +  — i  J  (31) 

d      l^p^     2  J 

The  heat  dQ  transferred  from  the  liquid  is  made  up 
of  two  parts:   (1)   the  convective  cooling  due  to  the 
temperature  difference  between  the  liquid  and  gas  and  (2) 
the  evaporative  cooling  due  to  the  latent  heat  supplied  to 
the  liquid  vaporized.   The  convective  cooling  is 


-dQ^  =  hA^N(T^  -  T  )dt  (32) 

2 
where  h  is  the  heat-transfer  coefficient.  A,  =  ttd   is  the 

d 

surface  area  of  a  drop,  and  dt  =  dx/V_  is  the  time  required 

Li 

to  traverse  dx.   Thus, 

6hm  (T   -  T, )dx 

-dZ   =  i — 2 i —  (33) 

c      Dp 

1^1 
The  evaporative  cooling  is 


-dQ^  =  L^dm^   +  L^di^  (34) 

g      g 
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The  change  in  enthalpy  of  the  liquid  over  the  pressure 
increment  dp  is 

dH  =  A,fc^dT^  +  3E  .  ^j  (35) 

Substituting  Eqs.  (31),  (33),  (34),  and  (35)  into  the 
steady-flow  energy  equation  dQ  -  dW  =  dH,  the  result  is 


6hin,  5Tdx 

— =;^ — — —  -  L  dm   -  L,dm,   =  m^CTdT,  (36) 

Dp,  V,      aa     h     Id  111 

11         g        g 

where  6T  =  T   -  T^ . 
g    L 

Writing  for  a  finite  interval,  the  final  form  of  the 
drop  heat-transfer  equation  is 

1 


At,  =  ■'- 


^   ^1 

m 


<n-  f  m  A        Am         Am, 
6h6  TAX         a  b 

m      ^ g  _ 


-  L    -r-^  -  L, 


Dpi  Vi      a   ihn      b   m, 
"^^m  ^m     ^   ^m      "^   ^m  j 


(37) 


Equations  (1) ,  (17) ,  (23)  ,  (28)  and  (37)  are  the  five 
equations  that  must  be  solved  simultaneously  to  obtain  the 
values  of  the  five  dependent  variables  D,  T  ,  T  ,  V  ,  and 
V-  as  a  function  of  the  independent  variable  p.   To  carry 

Li 

out  the  solution  all  quantities  in  the  equations  must  be 
expressed  in  terms  of  these  six  variables. 

The  preceding  equations  form  the  basis  for  the 
mathematical  model  which  is  used  to  predict,  based  on 
inlet  conditions,  the  exit  velocity,  and  temperature  of  the 
mixture.   These  equations  also  form  the  basis  for  the  model 
which  provides  the  optimum  nozzle  shape  given  a  set  of  inlet 
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conditions.   Some  additional  relationships  are,  however, 
required.   These  are: 

1.  Phase  properties — to  establish  the  mass  ratio, 
mass  flow  rate  ratio  of  gas  to  liquid,  and  the 
thermal  conductivity  of  the  mixture. 

2.  Liquid  drop  drag  coefficients. 

3.  The  liquid  drop  heat  transfer  coefficients. 

4.  Boundary  layer  momentum  thickness  and  displace- 
ment thickness. 

5.  Skin  friction  coefficient. 

These  five  additional  relationships  are  developed  in  detail 
in  Re f.  1. 

The  preceding  "basic"  model  in  conjunction  with  the 
additional  equations  form  the  core  for  four  computer  programs 
which  provide  real  and  isentropic  flow  output  parameters, 
and  an  optimized  nozzle  contour  given  a  set  of  inlet 
conditions. 

Reference  1  reported  that  the  optimum  nozzle  profile 
has  an  elongated  throat  region,  that  the  required  shape  is 
insensitive  to  drop  diameter  and  nozzle  length,  and  that  the 
exit  velocity  is  insensitive  to  departures  from  the  optimum 
shape.   For  given  fluids  and  pressures,  the  exit  velocity 
was  found  to  be  insensitive  to  nozzle  length  and  flow  rate 
beyond  certain  minimum  values.   This  insensitivity  of 
nozzle  exit  velocity  to  nozzle  shape  and  size  results  from 
the  inverse  dependence  of  drop  size  on  slip  velocity,  which 
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apparently  acts  to  adjust  the  slip  to  a  value  that  holds  the 
exit  velocity  in  the  range  75  to  90%  of  isentropic.   Thus, 
it  is  difficult  to  design  either  a  very  good  or  a  very  poor 
two-phase  nozzle. 
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V.   EXPERIMENTAL  APPARATUS 

The  experimental  apparatus  was  designed  to  provide  a 
facility  to  study  the  performance  of  a  Bi-phase  nozzle. 
The  basic  requirements  were  to  allow  overall  nozzle  per- 
formance evaluation  as  well  as  optical  diagnostic  probing 
of  the  Bi-phase  flow-field. 

The  experimental  system  can  be  conveniently  grouped  into 
five  subsystems.   These  are:   the  nozzle  proper ,  nozzle 
assembly,  air  supply  system,  liquid  injection  system,  and 
nozzle  exhaust  plane  instrumentation.   Each  subsystem  is 
described  in  the  following  sections.   Figure  8  is  an  overall 
system  drawing. 

A.   NOZZLE 

The  nozzle  was  designed  based  on  recommendations  given 

in  Refs.  1  and  5.   It  is  convergent-divergent,  12  in.  long 

2  2 

with  a  throat  area  of  0.25  in.  .   Inlet  area  is  2.0  m.  . 

The  throat  is  located  4  in.  from  theinlet  (1/3  of  the  over- 
all length) .   The  nozzle  is  constructed  by  "sandwiching" 
two  h   inch  thick  machined  aluminum  nozzle  profile  plates 
between  H   inch  plexiglass  plates  (Figs.  9  and  10) .   The 
aluminum  nozzle  plates  are  located  at  the  end  of  a  30  inch 
long  "test"  section  and  are  easily  replaceable.   The  test 
section  provides  the  necessary  space  for  pressure  taps, 
liquid  injection  tubes,  and  other  instrumentation  and 
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Fig.  8.   Overall  System  Drawing 
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NOZZLE  ASSEMBLY 


Inlet  Flange 
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.Plexiglas 


Aluminum  Plate 
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Fig.  9.   Nozzle  Assembly  Drawing. 
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Fig.  10,   Nozzle  Assembly  Photograph, 
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configuration  options.   The  ability  to  replace  the  nozzle 
profile  plates  allows  for  modification  of  the  nozzle 
configuration . 

B.   AIR  SYSTEM  (Fig.  11) 

3 
Two  117  ft.   tanks  supplied  with  compressed  air  from  an 

Ingersol  Rand  three  cylinder  air  compressor  provide  the  air 

storage  volume  required  to  support  nozzle  operation.   The 

nozzle  is  supplied  with  air  at  variable  pressure  via  3"  ID 

piping.   The  air  supply  piping  contains  two  pressure  control 

valves,  a  solenoid  actuated  nitrogen  operated  3  in.  ball 

valve,  and  a  flow  measurement  orifice  plate.   The  orifice 

flow  coefficient  was  determined  from  Ref.  6,  Table  4.   The 

air  mass  flow  rate  can  be  determined  from  the  equation 


W,  =  359  K  d  Fa  Y  /ETy 
h  w' 


where 


W,  -  air  mass  flow  rate  in  LBm/HR 

K   -  orifice  flow  coefficient  Table  4  [Ref.  6] 

d  -  orifice  diameter 

F   -  thermal  expansion  factor  Fig.  38  [Ref.  6] 
a 

y   -  expansion  factor.  Fig.  40b  [Ref.  6] 

h   -  differential  pressure  in  H^O  at  68  F 
w  ^2 

Y   -  Specific  weight  of  the  flowing  fluid  at  the  inlet 

3 

side  of  the  primary  element  (LB/FT  ) 
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Fig.    11.      Gas   Flow  Drawing. 
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and 


h  (y  -  Y  ) 
,      m  'm   'o 
h   = 


w     62.317 

where 

h   -  inches  of  differential  of  inches  of  manometer  fluid 
m 

3 

Y  -  specific  weight  of  manometer  fluid  LB/FT 

Y  -  specific  weight  of  fluid  separating  manometer  from 

flowing  fluid 
h   -  differential  head,  in.  of  H^O  at  68°F 

W  A 

And 

Y  =  Y^ 

Where 

Y  =  specific  weight  of  dry  air  at  actual  temperature, 
a 

For  an  example  calculation,  see  Appendix  2,    example  5  from 

Re  f .  6  . 

C.   LIQUID  INJECTION  SYSTEM  (Fig.  12) 

The  liquid  injection  system  consists  of  a  pressure 
vessel  to  contain  the  water  which  is  pressurized  from  a 
nitrogen  tank.   The  pressure  vessel  supplies  pressurized 
water  to  the  injection  tube  via  a  solenoid  stop  valve,  a 
flow  control  valve  and  a  flowmeter.   The  flowmeter  was 
calibrated  using  an  accurate  scale  and  a  stop  watch. 
The  calibration  data  and  chart  are  recorded  in  Fig.  13. 
The  liquid  injector  is  a  0.25  in.  brass  tube  inserted  in 
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Fig.  12.   Liquid  Flow  Drawing. 
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Fig.  13,   Liquid  Flowmeter  Calibration  Curve. 
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the  3  in.  ID  air  supply  pipe  just  upstream  of  the  flange 
connection  to  the  test  section.   The  injector  tube  is 
drilled  with  5  holes  facing  the  test  section  entrance. 
It  has  been  proposed  that  more  efficient  liquid  atomization 
could  improve  nozzle  efficiency;  therefore,  the  drilled 
holes  were  made  as  small  as  possible  consistent  with 
achieving  a  significant  liquid  mass  flow  rate. 

D.   NOZZLE  EXHAUST  INSTRUMENTATION 

A  cylindrical  drag  body,  which  is  attached  to  a  rectan- 
gular cantilever  beam,  is  placed  in  the  nozzle  exhaust 
flow  stream.   The  drag  body  unit  is  attached  to  a  bracket 
which  is  attached  to  a  rectangular  transversing  plate 
(Figs.  14  and  15) .   The  transversing  plate  is  screw  driven 
to  move  in  a  transverse  direction  across  the  nozzle  exhaust. 
This  movement  causes  the  cylinder  to  move  through  the 
entire  nozzle  exhaust  flow  stream.   By  placing  strain  gages 
on  the  drag  body  unit,  the  cantilever  beam  deflection  can 
be  determined  and  related  to  the  force  on  the  cylindrical 
portion  of  the  drag  body.   That  force,  in  turn,  can  be 
related  to  the  local  velocity  of  the  flow  stream  and  these 
velocities  plotted  to  give  a  velocity  distribution  at  the 
nozzle  exhaust.   A  similar  plot  can  be  constructed  for  the 
momentum  distribution 

„2 
Drag  (LB^)  =  C^^  A  ^ 
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Fig.    14.      Instrumentation   Drawing 
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Fig.  15.   Instrumentation  Photograph, 
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where 

C_.  -  Drag  coefficient  of  the  cylinder 

-  (Smooth  cylinder  1/8  in.  in  diameter) 
A   -  Projected  area  of  the  cylinder  1/2"  X  1/8"  =  1/16  in.^ 
p   -  Density  of  the  mixture 
U   -  Local  velocity  of  the  exhaust  flow  stream 

Solving  for  U: 


'2   X  Drag  force 

"=V — ^ — 


Figures  16  and  17  are  calibration  data  for  the  drag  cylinder. 
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Fig.    16.      Drag  Cylinder  Calibration   Curve. 
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Fig.  17.   Drag  Cylinder  Calibration  Curve 
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VI.   CONCLUSIONS 

Successful  Bi-phase  nozzle  operation  was  achieved  and 
sustained  under  controlled  conditions.   Liquid  injection, 
nozzle,  and  instriomentation  system  all  operate  satisfactori- 
ly.  The  liquid  breakup  and  acceleration  in  the  nozzle 
passage  was  observed. 

With  the  nozzle  pressures  employed,  i.e.,  20,  30,  40, 
and  50  psig  an  instability  was  observed  in  the  downstream 
portion  of  the  nozzle  passage.   It  was  apparent  that,  for 
the  nozzle  geometry  installed,  there  was  a  transition  from 
supersonic  flow  in  the  throat  region  to  subsonic  conditions 
in  the  diverging  downstream  passage.   It  is  apparent  that 
at  this  point  the  passage  converts  from  a  supersonic  nozzle 
to  a  subsonic  diffuser.   The  transition  point  appeared  to 
be  unsteady,  initiating  alternatively  from  the  sides  of  the 
nozzle.   The  position  at  which  transition  takes  place  could 
be  controlled  at  will.   As  the  liquid  mass  ratio  and/or 
upstream  pressure  were  increased,  the  passage  tended  to 
perform  as  a  supersonic  nozzle,  i.e.,  the  supersonic- 
subsonic  transition  was  driven  further  downstream  towards 
the  exit.   This  is  as  one  would  expect.   As  the  liquid 
mass  ratio  is  increased,  the  sonic  velocity  decreases  and 
the  passage  tends  to  a  supersonic  nozzle.   Also,  as  the 
pressure  ratio  increases,  the  shock  or  transition  to  subsonic 
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conditions  is  driven  downstream.   An  exit  plane  survey 
(Fig.  18)  using  the  instrumented  drag  cylinder  confirms 
the  existance  of  a  transition  point  somewhere  in  the  diverg- 
ing portion  of  the  flow  passage.   Thus,  one  is  led  to 
conclude  that  the  role  of  change  of  flow  passage  area  is 
too  high  for  the  test  pressure  ratios.   From  Fig.  18  it 
is  also  evident  that  the  flow  tended  to  oscillate  from  one 
side  of  the  nozzle  to  the  other.   This  was  confirmed  by 
visual  observations  as  well  as  audible  pulsations. 

It  would  be  desirable  to  develop  a  degree  of  flexibility 
in  changing  or  altering  the  nozzle  profiles,  particularly 
the  diverging  passage.   Originally  it  was  envisioned  to 
replace  the  nozzle  inserts;  however,  it  may  be  more  expedi- 
tious to  provide  an  infinitely  and  continuously  variable 
system.   This  may  be  achieved  by  a  flexure  system  controlled 
by  jack  screws  in  the  exit  plane.   The  second  desirable 
aspect  would  be  the  capability  to  diagnostically  probe  the 
flow  field  in  the  passage  proper.   A  laser  anemometer  may 
provide  the  capability;  however,  no  information  exists  on 
the  use  of  LDA  in  Bi-phase  flow. 
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